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Previous numerical investigations of an one-dimensional DNA model with an extended 
modified coupling constant by transcripting enzyme are integrated to longer time and 
demonstrated explicitly the trapping of breathers by DNA chains with realistic param- 
eters obtained from experiments. Furthermore, collective coordinate method is used to 
explain a previously observed numerical evidence that breathers placed far from defects 
are difficult to trap, and the motional effect of RNA-polymerase is investigated. 



1. Introduction 

The editors of Science, Culotta & Koshland [1994], announced in the last issue 
of 1994 that 'the molecule of 1994' is DNA. One can still find many papers in the 
physical literature about DNA recently. However, most physicists are concerned 
about its statistical properties by analyzing some database of DNA sequences, for 
instance Azbel' [1995] found DNA sequences have no long-range correlations, instead 
of dynamical properties analyzed below. 

An artistic plot of DNA which appeared on the cover page of that issue of 
Science and most textbooks of biochemistry is similar to Fig. |l|. As one can see, 
it is a double-helix with hydrogen bonds connecting those two strands. Of course, 
this is a simplified picture. However, Peyrard & Bishop [1989] untwisted this ladder 
like DNA, and wrote down its equation of motion. 

Biologists told us that transcription of RNA is the first step in a chain of events 
leading to expression of the genetic information encoded in double-stranded DNA. 
During this process RNA-polymerase will attached to DNA. Erie et al. [1993] 
consider the role of RNA-polymerase in transcription is to synthesize, under the 
DNA template, the nascent RNA chain with high fidelity and at reasonable rates. 
However, for physicists the enzyme may play an additional role before transcription, 
i.e. to focus thermal energy so that the double-helix can be opened. 

Englander et al. [1980] have proposed the existence of solitons on DNA helices 
more than 15 years ago. Technical difficulties prevent direct observation of solitons 
on DNA. However, Urabe & Tominaga [1981] have Raman spectra showing mode 
softening which might be accounted by the existence of breather traveling. On the 
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other hand, since DNA is a kind of polymer and solitons have been recognized as 
important conformational excitations in polymers. Furthermore, numerical and an- 
alytical investigations of the Peyrard-Bishop model of DNA, which can be rationally 
deduced from a real DNA atomic structure, show the existence of breathers travel- 
ing on DNA chains. We are now more confident about the existence of breathers 
on DNA chains today. 




Figure 1.: An artistic plot of DNA. 



Forinash et al [1991,1994] considered the interaction of enzyme with DNA during 
transcription as an isolated impurity. In a previous paper, Ting & Peyrard [1996] 
considered numerically the possibility of breather trapping by a continuous defect 
and by a multiple scale analysis the necessary condition for trapping to be a reduced 
coupling constant within the defect area. The conclusion was loose about whether 
the breather will be trapped by a real DNA. Due to limited computing resources, 
numerical integrations were only carried out to about 5000 time steps. However, 
Dauxois & Peyrard [1993] indicated that DNA is a discrete object and the behaviour 
of a breather at this region of coupling is mainly governed by discrete effect. In 
that case a stable breather should be thin and move slowly. It is the purpose of this 
paper to extend those numerical integration to longer time to investigate the long 
time behaviour of that kind of breather motion. The second purpose of this paper 
is to extend previous collective coordinate method to show the difficulty of trapping 
of a far away breather and to consider the effect of moving enzyme, because Erie et 
al. [1993] told us that RNA-polymerase is actually moving. 

2. The Peyrard-Bishop DNA Model 

The Peyrard-Bishop model considers a harmonic coupling for neighboring sites 
and within the same strand, while an anharmonic Morse coupling between strands. 
A transformation into normal coordinate is made and the longitudinal mode thrown 
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away, because it is too small. The equations obtained, in dimensionless form, read 

~ k n+1 (y n+1 - y n ) + k n (y n - y n _ x ) - 2e y "(e Vn - 1) = , (1) 

in which y n is the relative displacement between strands and k n is the coupling 
constant between y n and y n -i- Because we have untwisted the double-helix, there 
is no difference between A-, B- or Z- types of DNA. Furthermore, the present model 
does not taken specific DNA sequence into account. Ansari et al. [1995] told 
us, protein contact DNA at certain sites and bend DNA towards itself as shown 
in Fig. ||. This will results in a reduced coupling constant for the outside strand. 
Presumably, enzyme can function as a kind of lens to focus some thermo fluctuation 
waves present in normal physiological condition, traveling along the strand. 



Figure 2.: An artistic plot of a DNA with enzyme, showing the enzyme contact 
DNA at multiple sites and bend DNA towards itself. 



2.1. Long-time numerical integrations 

Previous numerical integrations can be summarized as the following: Taller 
breathers, hence thinner, can be trapped by a defect with reduced coupling con- 
stant, while shorter breathers are broader and result in only concentrated energy 
and pass through. Therefore trapping can be controlled by breather's amplitude. 
Furthermore, these two points indicate breathers move collectively as rigid bodies 
just like bike wheels. 

One may wonder whether the results obtained in the previous paper are appli- 
cable to more realistic cases. In real DNA one has k n rs 0.13 to 0.15. Two pictures 
for trapping by realistic DNA are shown in Fig ||. These breathers are so weak, 



4 DNA Transcription Mechanism 



i.e. move so slowly, so that they are vulnerable to the presence of defect or other 
breathers. In Fig. |^(a) the defect pattern is reversed in comparison with Fig. ||(a). 
But we found quite different picture for the breather motion and they were attracted 
by the region of lowered coupling constant. Fig. ^(b) shows two breathers, the same 
as in Fig. 0(a) except for their separation, attracting each other. 



0.2 
0.0 



1x10 



26 300 
site number 

5 



fx/0 



6x7 



4x10 



4 



2x10 




1 



00000 



90000 



^0000 



45000 



35000 



20000 



T 

I \ 10000 



1 45 



300 



1 150 300 



o 



Julian J.-L. Ting 5 



0.2 
0.0 



1x10 



26 300 
site number 

5 



fx/O 



6x7 



4x10 



2x10 



1x10 



Sx/0 




i 



Li 



4x70 



4 




-2x/0 



/ 45 



300 



I 100000 



75000 



60000 



K5000 



35000 



/oooo 



/50 300 



Figure 3.: Half-tone plots of the energy distribution of breather evolutions, left, and 
the corresponding amplitudes, right, from direct numerical integration of Eq. (|l|). In 
each figure the top insert shows the variation of the coupling constant used for the 
calculation, while the right inserts show snap-shots of the breather energy or am- 
plitude distributions. Defects positions are shown in the plots axis. The breathers 
start from X(0) = 130,240, K n = 0.15 outside the defect, the amplitude of sinu- 
soidal defect is 0.05 and (a) u c = 0.14, 0.13, u e = -0.05, -0.06. (b) u c = 0.13, 0.11, 
u P = -0.04,-0.03. 
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Figure 4.: Both figures are the same as Fig. ||(a), except (a) the defect pattern is 
reversed, (b) the initial breather position are closer. 



2.2. Collective coordinate method for moving defect 

Following Rcmoissenet [1986], a slightly different derivation from the one used 
before is used below to consider the effect of moving enzyme. Firstly, in the con- 
tinuum limit with a Taylor expansion in the potential term which assumes small 
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amplitude excitation Eq. ([TJ) became, 

Because of the small oscillation assumption, we have y ks e<j) and 

<j> = F^ 19 + F*e- W + e(F + F 2 e 2ld + F.*e- 2W ) + 0{e 2 ) , (3) 

in which F\ = F±{xi,t±) and F% = F2{x±,ti) with x\ = ex etc.. This expansion is 
equivalent to 

</> = F Q + eF 1 +e 2 F 2 + 0{e 3 ) (4) 

used before. Furthermore, for low frequency breathers we considered, we have 
6 w —Lit, and 

n d 

+ dx~ 1 € + 0{eh (5) 

— « — e (7) 

dx dxi 

for the smooth variation of the defect. Collecting terms of equal powers in e luJt , 



d 


d 


dx 


dx 


d 2 


d 2 


dx 2 


dx 2 



We also assumed 



F - 3F,F* = , 

* 2rj?F-LF*fiM 7 2 2 w * w 2^ h d 2 F l OF, 8 2 F l dk8F l 

3eauj {FoFi + Fi Fa) - r a u Fi Fi + fce ^ + 2,w - e w + & ^ = ' 

F 2 = --F 2 . 
2 1 

From the above three equations we obtained a perturbed Nonlinear Schrodinger 
equation (NLS) at order e 2 : 



2i m + 9k 1 8F^ +k ^ Q 

ot 2 dxi dxi ox{ 



We can rescale it into a standard form, 



1 1 d ~ 

iu t - + + u|u| 2 + -^r(fcu £ ) = , (9) 



and obtain the corresponding Lagrangian density, 



A - l -{u*u t - uu* t ) - 1(1 + fc)|t*| 2 + ~l«| 4 ■ (10) 
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A collective coordinate ansatz, 

u(x, t)=r) sech^x - Q , (11) 

brings the full space integrated Lagrangian density into 

n 3 1 f + °° 

L = -2ri<l>t - 2C& + -V - Z 2 V - r / fcKfdx. (12) 



3 * ' 2 _ _ 
For a moving defect 

k = K[B(x-vt + l)-G(x-vt-l)], (13) 
we have the equation of motion for the collective coordinate variables 



2 



-~(S£l*(i vt) + SlTUl + vt)W - j(Tl + W , (14) 

£t = ~(SlTl-S 2 _T 2 _) V s -^(Sl-S'_)ev, (15) 

Ct = ^+|(T++T_)^, (16) 

»?t = 0, (17) 

in which 

7+ = - tanh^ut - 77Z - () , (18) 

T_ = tanh(r/wt + ??/ - C) , (19) 

5+ = sech(r/ui - 77/ - () , (20) 

5_ = sech^vi + rjl - C) • (21) 

We find there is nothing essential different from previous stationary defect lo- 
cated at origin, which can be obtained as a special case of present result at t = 0, 
except for the variation of <fi t and the redefinition of Eqs. (|l8|)-(pT|). We can still 
prove that a reduced coupling constant in the defect area is a necessary condition 
for a breather to be trapped. Since a necessary condition for trapping is Ct — 
more than twice, which, according to Eq. (|l6|), is equivalent to 

cosh 2 (?yi;< — () = 1 — cosh 2 (?^) — Ksinh(f?Z) cosh(^Z). (22) 

Because cosh 2 > 1 and 77 > 0, k has to be negative. Furthermore, we can interpret 
the above equation as: the greater the separation between the breather and the 
defect; the larger the left hand side, the smaller k, i.e. the greater |«|, should 
be in order to balance the left hand side to fulfill the trapping requirement. This 
point is consistent with our previous numerical experience that if we put the initial 
breathers farther from origin we require larger breather amplitude for the breather 
to be trapped. However, the smallest k is —1. 
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3. Conclusion 

We have demonstrated explicitly the possibility for trapping of realistic DNA. 
Furthermore, we considered the effect of moving defect. The result shows no special 
influence on trapping through velocity except through changing the breather-defect 
separation distance under our collective coordinate model. However, a recent mea- 
surement by Yin et al. [1995] on the RNA-polymerase moving force might be 
taken into account by asymmetrical potentials and stochastic resonance, similar to 
Julichcr & Prost [1995] for molecular motors. 

There are other interesting works on DNA models: For instance, Yakushevich 
[1989] proposed another DNA model. More models can be found in the references 
of that paper and Forinash's [1991]. Within the Peyrard-Bishop model, Hisakado & 
Wadati [1995] have studied its inhomogeneity effect by considering random mass and 
random coupling constant. It will also be interesting to take long range columbic 
interaction into account as Gaididei et al. [1995] considered. Slepyan et al. [1995] 
considered solitary waves traveling on a three dimensional helicoidal fiber, which 
might be extendible to a three dimensional DNA model. However, much work 
remain to be done in the future. 
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